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Abstract 

Differential and total breakup cross sections as well as vector and tensor 
analyzing powers for p-d scattering are studied for energies above the deuteron 
breakup threshold up to Ei a b = 28 MeV. The p-d scattering wave function 
is expanded in terms of the correlated hyperspherical harmonic basis and 
the elastic S'-matrix is obtained using the Kohn variational principle in its 
complex form. The effects of the Coulomb interaction, which are expected to 
be important in this energy range, have been rigorously taken into account. 
The Argonne AV18 interaction and the Urbana URIX three-nucleon potential 
have been used to perform a comparison to the available experimental data. 
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I. INTRODUCTION 



In ref. the authors recently presented an application of the Kohn Variational Principle 
(KVP) in its complex form to calculate the elastic observables in p-d scattering for ener- 
gies above the deuteron breakup threshold (DBT). The KVP was implemented to describe 
continuum states of three outgoing particles including the distortion due to the Coulomb 
interaction in the asymptotic region. Only two energies were considered, Ei a i, = 5 and 10 
MeV. The validity of the KVP for the elastic S'-matrix describing the 2 — > 2 process in p-d 
scattering for energies above the DBT has been extensively discussed in ref. [[|. 

In the present paper the analysis of the elastic p-d reaction is extended up to Ei a b = 28 
MeV, covering the region where Coulomb effects are expected to be important. The large 
amount of accurate experimental data allows for interesting comparisons. It should be noted 
that, at present, the analysis of the polarization data at energies above the DBT has been 
done mainly by comparing n-d calculations to p-d data |||4j]. Differential cross section and 
vector analyzing power data exist for both n-d and p-d scattering, allowing for an estimate 
of the Coulomb effects. Conversely, no n-d data are available for the deuteron analyzing 
powers zTii, T 20 , T 2i , T 22 . These quantities are evaluated from experiments using a polarized 
deuteron beam on unpolarized proton targets. The inverse experiment of unpolarized proton 
or neutron beams on a polarized deuteron target seems to be extremely difficult at low 
energies and has not yet been done. 

Experiments using charged particles are certainly easier to perform and show smaller 
error bars than those using a neutral beam. On the other hand, the theoretical description 
of collisions with more than one charged particle in the final state has represented a difficult 
problem for many years. In a recent work a complete solution of the reaction e~ + H — > 
H + + e~ + e~ has been obtained by Rescigno et al. || by transforming the Schrodinger 
equation using the so-called exterior complex scaling and making use of supercomputers 
to solve the associated equations numerically. This was the first complete solution of a 
three-body collision with all the charged particles moving away from each other in the 
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final state. Regarding the p-d reaction, different techniques have been applied so far. The 
Faddeev equations in momentum space have been adapted to take into account the long- 
range Coulomb interaction using the screening and renormalization approach ||. Recently, 
a detailed comparison between the solutions of the Faddeev equations in configuration space 
and the KVP has been performed, though restricted to energies below the DBT [J/J. In the 
present work we turn our attention to describing p-d elastic observables above the DBT. In 
this case the application of the KVP is feasible and the calculation of the elastic S'-matrix 
does not require large computational devices. 

The study of the three-nucleon (3N) continuum provides important information about 
the capability of modern NN potentials to describe the three-nucleon dynamics. At present, 
a few realistic NN potentials are available that reproduce a large set of two-nucleon (2N) data 
with x 2 ~ 1 (per datum). They are substantially equivalent in reproducing all the details 
of the NN scattering, but in the description of nuclear systems with A > 2 differences 
appear. In addition, the three-nucleon system is the simplest one in which three-nucleon 
force (3NF) effects can be studied. The first signal for the necessity of a 3NF comes from 
the underbinding of the triton when only NN forces are used. Widely used 3NF models are 
based on the exchange of two pions with an intermediate A excitation. In general these 
models include a certain numbers of parameters which are not precisely determined by 
theory, so some of them can be taken as free parameters in order to reproduce, for example, 
the triton or 3 He binding energy. As a consequence, other observables which scale with the 
three-nucleon binding energy improve as well. Examples are the bound state r.m.s radii and 
the zero energy total cross section in n-d and n- 3 H scattering. On the contrary, vector and 
tensor N-d analyzing powers do not present such a scaling. 

Accurate measurements of p-d observables below the DBT, have been reported re- 
cently |8|-|T0| . A comparison of the theoretical predictions to these data shows an underpre- 



diction of the deuteron vector analyzing power iTu by ~ 30% [|TTJ. A similar discrepancy 
had been observed earlier in the neutron analyzing power A y , a problem which is usually 



known as the A y puzzle ||12|| . As the energy increases the observed discrepancies in A y 



and zTn reduce and tend to disappear, though not completely, above 30 MeV 0. Accord- 
ingly, the study of these observables over the energy region considered here is important for 
understanding such a behavior. 

Accurate 3N and 4N scattering wave functions are necessary for calculating a number 
of nuclear reactions. The technique used in the present work is based on the expansion 
of the wave function in terms of Jastrow type Correlated Hyperspherical Harmonic (CHH) 
basis functions. When the correlation factor reduces to a pair correlation function the 
Pair Correlated Hyperspherical Harmonic (PHH) basis is obtained. The CHH and PHH 



bases have been used to calculate the bound states of the A = 3,4 nuclei flTjyij] , N-d 
scattering f[5,|l6|, and p- 3 He and n- 3 H scattering [17| at energies below the three-body 



fragmentation. Moreover, wave functions obtained through those expansions have recently 
been used to study the radiative capture p + d — > 3 He + 7 below the DBT |Tj| and the hep 
process, namely the weak capture p + 3 He — > 4 He + e + + v e at the Gamow peak |19| . These 
two reactions have considerable astrophysical relevance. The former is the second reaction 
in the pp solar chain and has a prominent role in the evolution of protostars whereas the hep 
process plays an important role in the solar neutrino problem. The calculation of the p-d 
wave functions above the DBT will provide the input for further studies of radiative capture 
and photo and electrodisintegration of 3 He. 

In the present paper we present the results obtained for the differential and total breakup 
cross sections, nucleon analyzing powers A y and deuteron analyzing powers iT u , T 20 , T 2i 
and T 2 2 for N-d scattering at different energies. The calculations have been done using the 
two- nucleon AV18 potential p0| with and without the three-nucleon URIX force pl |. The 



results are given at nine different energies in the range 5 MeV < Eiab < 28 MeV. It has 
to be noted that, disregarding small corrections, E cm ^ = \E N {\E d ), where E N (E d ) is the 
nucleon (deuteron) incident energy and in the following we define E[ a b = E^. The highest 
energy considered here is E d = 56 MeV (Ei a b = 28 MeV) at which the deuteron analyzing 
powers are available [021. Just above the DBT, deuteron vector and tensor analyzing powers 



are available at i^ & = 5 MeV p3 |. For E^ < 18 MeV differential cross sections, proton and 



deuteron analyzing powers have been measured at several energies ||24j| . In ref. |25| differential 
cross section as well as vector and tensor observables have been measured between 8.5 MeV 
< Eiab < 22.7 MeV, though data for T21 are missing at some energies. 

The paper is organized as follows. In Section II the Kohn variational principle is reviewed. 
In Section III the numerical solution of the related differential equations are compared to 
previous results. Cross sections and observables are compared to the data in Section IV, 
and the conclusions are given in the last section. 

II. THE KOHN VARIATIONAL PRINCIPLE ABOVE THE DEUTERON 

BREAKUP THRESHOLD 

In the literature several investigations regarding the validity of the KVP above the DBT 
can be found, starting with the works of Nuttall |26| and Merkuriev ||27|| , where the dis- 



cussion, however, was limited to the n-d reaction. The first extensive demonstration of the 
applicability of the principle to the p-d collision has been given in ref. 0. The main result 
derived in |2j is that the effect of the Coulomb interaction can be taken into account in 
such a way that the form of the principle remains unchanged when the energy goes from 
below to above the DBT. Below the DBT the collision matrix is unitary and the problem 
can be formulated in terms of the real reactance matrix (if-matrix) . Above the DBT the 
elastic part of the collision matrix is no longer unitary and the formulation in terms of the 
S- matrix, the complex form of the KVP, is convenient. Refering to ref. || for details, a brief 
description of the method is given below. The scattering wave function (w.f.) $ is written 
as sum of two terms: 

* = *c + ^a- (1) 

The first term, ^c, describes the system when the three-nucleons are close to each other. 
For large interparticle separations and energies below the DBT it goes to zero, whereas 
for higher energies it must reproduce a three outgoing particle state. It is written as a 
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sum of three Faddeev-like amplitudes corresponding to the three cyclic permutations of the 
particle indices 1, 2 , 3. Each amplitude ^(x,, y;), where x», y, are the Jacobi coordinates 
corresponding to the i-th permutation, has total angular momentum JJ Z and total isospin 
TT Z and is decomposed into channels using LS coupling, namely 

*c(x», y<) = E Mxi, yi)y a (Jk, i) (2) 

a=l 

y a (jk,i) = {[Y ia {mLM\A a [^j*\ Sa } JJz \€t l a ] TTz , (3) 

where x iy yi are the moduli of the Jacobi coordinates and y a is the angular-spin-isospin 
function for each channel. The maximum number of channels considered in the expansion 
is iV c . The two-dimensional amplitude <p a is expanded in terms of the PHH basis 



<fra(Xi,yi) = P 5/2 fa(Xi) 



E<(p) (2) ^t Ld 



(4) 



L K 

where the hyperspherical variables, the hyperradius p and the hyperangle fa, are defined 
by the relations Xi = p cos fa and yi = psin0j. The factor ^ P^ L (</)) is a hyperspherical 
polynomial and f a (xi) is a pair correlation function introduced to accelerate the convergence 
of the expansion. For small values of the interparticle distance f a (xi) is regulated by the 
NN interaction whereas for large separations the correlation function is chosen to satisfy 



f a ( Xi ) -> 1 0. 



The second term, ^/a, in the variational wave function of eq.fl]]) describes the asymptotic 
motion of a deuteron relative to the third nucleon. It can also be written as a sum of three 
amplitudes with the generic one having the form 

n x LSJ (x l ,y l )= £ w^nliy,) {{[YU^s^s*] Y^)} [t^] TTz , (5) 

l a =0,2 b JJz 

where Wi a (xi) is the deuteron w.f. component in the state l a = 0,2. In addition, = 
1, t J a = and L is the relative angular momentum of the deuteron and the incident nucleon. 
The superscript A indicates the regular (A = R) or the irregular (A = 7) solution. In the 
p — d (n — d) case, the functions TZ X are related to the regular or irregular Coulomb (spherical 
Bessel) functions. The functions Q x can be combined to form a general asymptotic state 
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Yi) = ^L5j(x l; Yi) + E ^flX'SV^ Y<) > ( 6 ) 



L'S' 



where 



fi Lj( x i> Yi) = Moo^L 5 j( x i> Yi) + «oi^£ 5 j(xi, Yi) , ( 7 ) 
fi L?j( x i> Yi) = Mio^L 5 j( x i> Yi) + M ll^LSj( x i> Yi) ■ ( 8 ) 

The matrix elements be selected according to the four different choices of the matrix 

C = iT-matrix, i^ _1 -matrix, S-matrix or T-matrix. A general three-nucleon scattering w.f. 
for an incident state with relative angular momentum L, spin S and total angular momentum 
J is 



Hsj= E [M*>yi) + n£sj(*i>yi)] > ( 9 ) 



i=l,3 

and its complex conjugate is ^isj- A variational estimate of the trial parameters in the 
w.f. *&lsj can be obtained by requiring, in accordance with the generalized KVP, that the 
functional 

['£!*] = J L% - ^(^lsjIH - EWhpj) , (10) 

be stationary. Below the DBT due to the unitarity of the S'-matrix, the four forms for the 
£-matrix are equivalent. However, it was shown that when the complex form of the principle 
is used, there is a considerable reduction of numerical instabilities |p8| . Applications of the 



complex KVP for N-d scattering (below the DBT) can be found in ref. |29j . Above the DBT 
it is convenient to formulate the variational principle in terms of the S'-matrix. Accordingly, 
we get the following functional: 

[ J si$] = J s s L s L : + i(nsj\H - E\n, s >j) ■ w 

The variation of the functional with respect to the hyperradial functions u^{p) leads to 
the following set of coupled equations (hereafter named SE1): 

E V&Wjz + B w^To + c ^' {p) + 1F E N ™^\ u «^ = D ^ ■ (12) 

ct'k f i > 



For each asymptotic state ( 2S+1 *>Lj two different inhomogeneous terms are constructed corre- 
sponding to the asymptotic Q^sj functions with A = 0, 1. Accordingly, two sets of solutions 
are obtained and combined to minimize the functional (|ll|) with respect to the S'-matrix ele- 
ments. This is the first order solution, the second order estimate of the S-matrix is obtained 
after replacing the first order solution in eq.(|ll]) fL5U29f . 

In order to solve the system SE1 appropriate boundary conditions must be specified 
for the hyperradial functions. For energies below the DBT they go to zero when p — > oo, 
whereas above the DBT energy they asymptotically describe the breakup configuration. The 
boundary conditions to be applied in this case have been discussed in refs. @,[3(J and are 
briefly illustrated below. To simplify the notation let us label the basis elements with the 
index p = [a, if], and introduce the following completely antisymmetric correlated spin- 
isospin-hyperspherical basis elements 

V,(p,Q) = j^Uxi) i2) PK a ' La (&)y a (jk,*) , (13) 

8=1 

which depend on p through the correlation factor and form a non-orthogonal basis. In terms 
of the Vfj,(p, £1) the internal part is written as 

y c = p- 5/2 Y.Mp)'P»(p,n) , (i4) 

with N m the total number of basis functions considered. The "uncorrelated" basis elements 
V®(Q) are obtained from eq. (|H|) by setting all the correlation functions f a (xi) = 1. It is 
important to note that the elements V®(£l) do not form an orthogonal basis, as has been 
discussed in ref. |JT|] where the standard hyperspherical harmonic basis (HH) has been used 
to calculate the three-nucleon bound state. Those basis elements having the same grand- 
angular quantum number = £ a +L a +2K, the same A a and S a , but belonging to different 
channels, are not orthogonal to each others. Moreover, some of them are linearly dependent. 



In Ref. [^1|] such states have been identified and removed from the expansion used to describe 
the triton bound state. 

In the present case, the basis elements V^{p — > oo,Q) reduce to the uncorrelated ones 
V®(Q) in the asymptotic region since f a {x) — > 1 for large interparticle distances. Therefore, 
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it appears useful to combine the correlated basis ([13]) in order to define a new basis with the 
property of being orthonormal when p — ► oo. This can be readily accomplished by noting 
that the matrix elements of the norm N behave as 

, jV (3) 

iV(p) = / dn v,( P , n)t^(p, n) - iv$ + -f. + o(i/ P 5 ) , for P -> oo , (15) 

where, in particular, 

= f<mr°(nyv°,(n). (ie) 

Let us define a matrix {7 such that the matrix Z7* TV^ ^ U — J\f is diagonal with diagonal 
elements A/" M either 1 or 0. The values A/" M = correspond to states "P°(fi) that depend 
linearly on others. New uncorrelated and correlated bases are defined as: 

Qj(o) = £ u^v°,(n) , q,( p , n) = £ tWVfa fi ) . ( 17 ) 
^'=i /i'=i 

The basis functions Q^p, O) are still not orthogonal for any finite values of p. When p — > oo, 
the elements Q M (p, fi) — > Q„(0). Due to the fact that some of the uncorrelated elements 
V®(Q) are linearly dependent, some elements Q°(fi) are identically zero. Therefore, some 
correlated elements have the property: Q^(p,Q) — > as p —>■ oo. In the following we arrange 
the new basis in such a way that for values of the index p < N m the eigenvalues of the norm 
are Af^ = 1 and for N m + 1 < p < N m they are Af^ = 0. 
In terms of the new basis, the internal part is simply 

*c = p - 5/2 Y,Mp)Q»(p^) , (is) 

where the old set of hyperradial functions is related to the new set through the transformation 
= Ylu' U^i coy . The variation of the functional ( |10D with respect to the new hyperradial 
functions w^p), which are now the unknown quantities entering into the description of 
the internal part of the w.f. \I/cs leads to a set of inhomogeneous second order differential 
equations formally equal to SE1, and hereafter called SE2, in which each matrix X = 
A,B,C,N of eq.(jr2|) is substituted by X = U f XU and the inhomogeneous term D x by 
TJ A = U t D x . 
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For p — > oo, neglecting terms going to zero faster than p 2 , the asymptotic expression of 
SE2 reduces to the form 

^{-l (£ - ^ + ^) + ^ + °c?)Kw ^ ° • < 19 > 

where = h 2 Q 2 /M^, JC^ = + 3/2 and the matrix x is defined as 

= / dn Qj(n)t x Qj,(n) . (20) 

The dimensionless operator x originates from the Coulomb interaction as 

„ _ Mjy s e 2 1 + r j|Z 1 + r k)Z 

X 2^ 2 Q^cos0 i 2 2' [ } 

It should be noticed that Xw*' 

= if p,p' > iV m . 

In practice, the functions ^(p) are chosen to be regular at the origin, i.e. cu M (0) = and, 
in accordance with the equations to be satisfied for p — > oo, to have the following behavior 
(p < N m ) 

Nm 

- - £ (e-^ ln2Q " , V e lQp , (22) 
^'=1 



where £y are unknown coefficients. This form corresponds to the asymptotic behavior of 
three outgoing particles interacting through the Coulomb potential [|3~2" 1. In the case of n — d 



scattering (x = 0) the outgoing solutions evolve as outgoing Hankel functions H^ l \Qp) 

For values of the index p > N m the eigenvalues of the norm are A/" M = and the leading 
terms in eq.([0J) vanish. So, the asymptotic behavior of these functions is governed 
by the next order terms. A lengthly analysis of the 1/p 3 and 1/p 4 terms for each matrix 
X = A,B,C,N shows that these functions behave as e*^ p ~ E,tIn2< ^ where the quantities 
Q'^, £ M are related to the asymptotic expansion of the matrices A,B,C, N. This asymptotic 
behavior has been obtained neglecting all couplings between the p-th equation (p > N m ) 
and all the others. If couplings up to 1/p 4 are taken into account the quantities Q', E become 
matrices and we have (p > N m ) 
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MP) - " E eW'^ sln2 ^ , cy , (23) 
# -i mm 

where the cy are unknown coefficients. Previously we have shown that, for p > N m , the 
elements — > as p — > oo. The specific form of the (complex) matrix £ is such that in 
all cases oj^Q^ — ► as p — > oo. Accordingly, the states with p > iV m do not contribute to 
the outgoing flux. 



In ref. J3(J the set of equations SE2 has been solved numerically by choosing a grid 
of values for the hyperradius from the origin up to a certain value po- The differential 
operators have been substituted by finite differences in such a way that SE2 reduces to a set 
of linear equations that can be solved by standard numerical methods. In order to completely 
determine the problem, boundary conditions must be imposed at p = po- To accomplish 
this, eq.(|T§D has been solved for p > p taking into account coupling terms up to p~ 4 by 
an expansion of the functions in powers of 1/p and verifying the outgoing boundary 



conditions of eqs.(^J,^). Then, the continuity of the solutions and their first derivatives 
has been imposed at the matching radius po- The value of po is not important provided 
that the asymptotic expression of SE2 is already reached. This condition is well verified 
for values of the matching radius po ~ 80 — 100 fm. However, the functions ^(p) show 
an oscillatory behavior outside the range of the potential, typically for hyperradial values 
p > 30 fm. Therefore a large number of grid points were necessary to obtain stable solutions. 
Thus, in ref. [[!(| the calculation of N-d scattering states above the DBT was restricted to a 
simplified interparticle potential, namely an s-wave interaction. In such a case the number 
of coupled equations to be considered was sufficiently small. When realistic NN interactions 
are considered the number of coupled equations to take into account increases considerably. 
As a consequence, the dimension of the matrices after the reduction of derivatives to finite 
differences can be quite large. In order to keep the dimension of the matrices low, an 
alternative method of solution in the region p < p is to expand the hyperradial functions 
in terms of Laguerre polynomials [fij plus an auxiliary function 
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M 

MP) = P 5/2 E W exp(-f ) + < +1 ^(p) , (24) 

m=0 

where z = 7p and 7 is a nonlinear parameter. The linear parameters A™ (m = 0, M+ 1) 
are determined by the variational procedure. The functions defined above are matched to 
the outgoing solutions of eq. ( |I9|) at p = p . 



The inclusion of the auxiliary functions uJ^(p) defined in eq.(|24"D is useful for reproducing 
the oscillatory behavior shown by the hyperradial functions for p > 30 fm. Otherwise a 
rather large number M of polynomials should be included in the expansion. A convenient 
choice is to take them as the solutions of a one dimensional differential equation correspond- 
ing to the p-th equation of SE2: 

AMj^ 2 +BM- + CM + Q 2 WM]WM = I>>) . (25) 

The functions ID^ are chosen to be regular at the origin and they are matched to the solu- 
tions of eq.(^) which have been obtained through an expansion in inverse powers of p as has 
been previously discussed. For p > N m the matching at p has been done disregarding the 
couplings between the different equations in the region p > p 0l i.e. oJ^p) — > e l<yQ '^ p ~ T '' iln ' 2 ^ p \ 
As stated before, these states do not contribute to the outgoing flux and their importance 
in the construction of the scattering state diminishes very rapidly for large values of p. The 
approximation introduced for p > p in the application of the boundary condition to the 
states with p > N m has been checked by increasing the value of the matching radius. In 
the cases considered here the solutions obtained for the S'-matrix show a complete stability 
for values of the matching radius po > 100 fm. 

Let us define \p, m > to be a correlated totally antisymmetric element of the expansion 
basis. Here p indicates the correlated HH state Q^p, fi) and m = 1,....,M indicates the 
Laguerre polynomial L$(z) or, for m = M + 1, the auxiliary function aJ M . In terms of these 
basis elements the internal part of the wave function is 

* c = 5»,m> . (26) 

fi,m 
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The variation of the functional [ J Slf j ] with respect to the linear parameters leads to the 
following set of linear equations 

£ A$ <^m\H- E\im', mf >= L% m , (27) 

where the inhomogeneous term is 

Km = £<//, H# " E\Q X LSJ (^, y,) > . (28) 

3 

The first order solution of the S'-matrix is obtained solving the algebraic equations 

E J s a s x^,: = Y *g , (29) 

L"S" 

with the coefficients X and Y defined to be 

yffi =< n° LS j + Ksj\H - E\n° LSJ > , (30) 

where ^lsj * s constructed using the solution of eq.(^) with the corresponding inhomoge- 
neous term. The second order estimate [ J S££,'] is obtained replacing the first order solution 
in eq.(|TTJ). 

Finally, due to flux conservation the following condition has to be satisfied between the 
matrix elements of the elastic S-matrix and the coefficients of the outgoing breakup waves: 

El J sff'l 2 + £iM 2 = i- (si) 

S'L' ft 



The coefficients b^, which are defined in eq. ( 22 ) , are the linear parameters Aff +1 of eq.(^). 
The above relation allows the calculation of the total breakup cross section from the elastic 
S'-matrix elements, as has been recently discussed in ref. |34[ . 



III. NUMERICAL RESULTS 



In order to study the solution of eq.(|T2"D by means of the expansion given in eq. Q2~3]) , we 
have first calculated the phase-shift and inelasticity parameters for n-d and p-d scattering 
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using the spin-dependent s-wave potential of Malfliet and Tjon. The results are presented 
in Table I for two energy values, = 14.1 and 42.0 MeV. The calculations have been 
done using N a = 8 hyperspherical polynomials per channel, as in the case already studied 



in ref. []30| where the set of equations SE1 was solved using the finite difference technique. 
Moreover, since the potential is central, the phase-shifts 25+1 5l and inelasticities 2S+1 i]l do 
not depend on the total angular momentum J. Only the case L = has been considered and 
the results are given in Table I for increasing values of the number of Laguerre polynomials 



M. For the sake of comparison the results of ref. |3(J are reported as well as the benchmark 
results of ref. |35| obtained by solving the Faddeev equations in configuration space (Los 
Alamos group) and momentum space (Bochum group). We observe a very fast convergence 
with M and, in general, 16 to 20 polynomials are enough to obtain the phase-shift and 
mixing parameters with four digit accuracy. With the number of Laguerre polynomials that 
has been taken into account a very low dependence on the nonlinear parameter 7 has been 
observed. In fact the results reported here do not change for variations of the parameter in 
the range 1.5 fm -1 < 7 < 2.5 fm _1 . Moreover, the dimension of the matrices involved in 
the solution is one order of magnitude smaller than that used in ]30j . 

The case of realistic interactions has been considered in refs. [TJ|| where the AV18 in- 
teraction has been used to calculate p-d scattering at Ei a b = 5 and 10 MeV. In particular, 
in ref. [§] the convergence of the phase-shift and mixing parameters for the state J = l/2 + 
has been studied by increasing the number of angular-spin-isospin channels. The convention 
discussed in ref. |36| has been adopted in the parametrization of the S'-matrix in terms of 
phase-shift and mixing parameters. In order to illustrate the variation of these parameters 
with energy the doublet and quartet S, P and D phases, denoted as 2S+1 Lj, are reported 
in Fig.l as well as the mixing parameters 771/2+ , 773/2+) C1/2- and e 3 / 2 _. Both the real and 
imaginary parts are shown. It is interesting to notice that the splitting in the real part of the 
phases with equal spin S and angular momentum L but different J, increases with energy. 
Conversely, the imaginary parts of the phases, which are related to the inelasticity of a state 
with a given value of J, reveal a tiny splitting. After summing all the contributions, the 
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total breakup cross section can be obtained, as is discussed in the next section. 



IV. P-D CROSS SECTIONS 

The calculation of scattering observables using the present variational technique is based 
on the estimate of the elastic S- matrix for all states with J < J Max- Each observable is 
obtained from a trace operation after the evaluation of the transition matrix, following the 
formalism of Seyler |37| . The value of J Max has been chosen by requiring that partial waves 
with J > J Max give negligible contributions to all the observables considered. In the present 
work results for cross sections, vector and tensor analyzing powers up to = 28 MeV are 
presented, and correspondingly the value J Max = 19/2 has been found to be appropriate. 

Let us start with the analysis of the p-d cross sections. For p-d scattering the total 
breakup cross section accounts for all possible configurations in which all three particles 
are moving away from each other. Its expression can be given in terms of the elastic S- 



matrix [33] 



<P - d ) = ^\ E( 2J + " Sj&j}, (32) 

where k 2 = 2fiE cm /h 2 (fx is the nucleon-deuteron reduced mass) and Ij is the 3x3 identity 
matrix, except for J = 1/2 which is the 2x2 identity matrix. The quantity Sj is the elastic 
^-matrix for the state J. The sum runs over all possible values of J and parity (the sum over 
the two parities is implied). In principle the sum runs from J = to infinity, but there is a 
rapid convergence since each Sj matrix becomes closer to unitary as J increases. In Fig. 2 
the theoretical prediction for <Jb(p — d) is given together with the two sets of data available 



in the literature. The first data set corresponds to energies just above the DBT [38] whereas 



the second starts at 20 MeV ||39|| . The solid line is the AV18 prediction and is found to be in 
reasonable agreement with both sets of data. The inclusion of the URIX potential does not 
produce appreciable modifications and both results, with and without the inclusion of the 
3NF, nearly coincide. The low sensitivity to the 3NF can be understood by noticing that 
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the contribution to comes from a balance between the spin factor 2J + 1 and the quantity 
tr{Ij — SjSj} which can be considered as a measurement of the inelasticity of the state 
(divided by tr{Ij}). Above 5 MeV the state J = 3/2~ gives by far the main contribution to 
the observable |[34|| . The state J = l/2 + , which is appreciably modified by the 3NF, has the 
largest inelasticity, but due to a small spin factor it gives a contribution of the same order 
as other states that are much less "inelastic" and modified slightly by the 3NF. The final 
result after summing up all these contributions is that the small (but sizeable) effect of the 
3NF on J = 1/2+ has no impact in o&. 

Regarding the elastic p-d differential cross section, a huge amount of high quality data 
has been collected during the past years. Low energy measurements have been taken re- 
cently at TUNL at different energy values below Ei ab = 1 MeV |]ll3| , [4T||4"T|| . An analysis of 



the quality in the description of these data has been performed using the AV18 and the 
AV18+URIX interactions [f42| , |43"l| . It was shown that 3NF effects can be revealed through a 
X 2 analysis of the data. Essentially these effects are related to a correct description of the 
3 He binding energy. In fact, using the AV18+URIX interaction it is possible to describe 
the p-d differential cross section at E[ ab = 1, 2 and 3 MeV with a x 2 P er datum (xlf) 
close to one. This value increases significantly when the AV18 potential is considered alone. 
The agreement between the theoretical and experimental differential cross section worsens, 
though not dramatically, as the energy increases. For example, at Ei ab = 135 MeV a value of 
Xn = 16.9 (225.2) was recently obtained with (without) the inclusion of a 3NF [44]]. Again 
the inclusion of a 3NF reduces the x 2 P er datum considerably. 

The results obtained for the p-d differential cross section are given in Fig. 3 for nine values 
of the energy, Ei a b = 5, 7, 9, 10, 12, 14, 16, 18, 22.7, 28 MeV. For each energy three curves are 
shown corresponding to calculations using the AV18 potential (solid line), the AV18+URIX 
potential (dotted line), and calculations for n-d scattering using the AV18 potential (dashed 
line). The theoretical predictions are compared to the experimental data of refs. p4] , p5| , |22| , 
with the exception of the calculations at 16 MeV which are compared to data obtained at 
a slightly different energy (16.5 MeV). The analysis of the results at the different energies 
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shows that 3NF effects are small in this energy range and the AV18 and AV18+URIX curves 
practically overlap each other. A more quantitative analysis at Ei ab = 18 MeV gives x% = 11 
using the AV18 interaction and nearly the same value for AV18+URIX. Coulomb effects are 
mainly observed at forward and backward angles whilst they are strongly reduced at the 
minimum. Tiny Coulomb effects at the minimum are confirmed by comparing n-d to p-d 
data, as can be seen in ref. ||. 

As far as the agreement between theory and experiment is concerned, the situation for 
the differential cross section above the DBT is different when compared to what has been 
observed below the DBT. As mentioned before, at very low energies the differential cross 
section can be described with x% ~ 1 using AV18+URIX, while when the AV18 is used alone 
a substantially worse result, x% > 10; * s obtained. In fact, the AV18 curve remains above 
the data points all over the angular distribution. As the energy increases the tendency for 
the AV18 curve is to go below the data at the minimum. This problem is appreciable already 
at 28 MeV, as can be seen in the last panel of Fig.3. Around 20 MeV the AV18+URIX 
curve starts to rise above the AV18 curve and closer to the data. This effect is clearly 
shown in ref. Q where Faddeev calculations using several NN and 3NF interactions have 
been compared to the data at 3, 65, 135 and 190 MeV. In the energy range analyzed here 
we observe that there is one energy, around 18 MeV, where the AV18 and AV18+URIX 
curves mostly overlap. In order to analyze further this behavior, in Table II, the values at 
the minimum of the p-d cross section calculated with AV18 and AV18+URIX are compared 
to the data. The corresponding values of the AV18 n-d cross section are also given in order 
to have a quantitative idea of the size of the Coulomb effects. 

V. POLARIZATION OBSERVABLES 

The vector and tensor analyzing powers are examples of polarization observables. There 
is a large amount of p-d and d-p data for the vector analyzing powers A y and iTn as well 
as for the tensor analyzing powers T 20 , T 2 \, T 22 . The study of these observables is important 
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because they are sensible to the non-central terms of the nuclear interaction. These terms 
are responsible for small components in the wave function which in general are less known. 
Therefore, the accuracy shown by the modern interactions when reproducing the vector 
and tensor analyzing powers in the three-nucleon system gives important information about 
parts of the nuclear interaction not completely under control. As is well known, in the low 
energy region the vector analyzing powers are heavily underpredicted by all modern NN 
interactions and the origin of this discrepancy is not yet completely understood. Possible 
ways for solving this puzzle have recently been investigated, based on the inclusion of new 
terms in the three-nucleon potential |^5| , ^6| or on a new NN potential obtained from chiral 
perturbation theory |47[] . These studies represent only a first step in the understanding of 



the puzzle and further investigations and refinements of the models are needed. A similar 
underprediction of the proton analyzing power A y has been found in calculations on p- 3 He 



scattering, as was recently pointed out Therefore, a solution to the puzzle should 

concern both the 3N and 4N systems. 

In the present paper, we will discuss the quality of the description of the vector and 
tensor polarization observables achieved by the AV18 and the AV18+URIX interactions in 
p-d scattering up to 28 MeV. In Fig.4 the results for A y are given for the same nine energy 
values given in Fig.3. The three curves correspond to the p-d A y calculated using AV18 
(solid line) and AV18+URIX (dotted line), and the n-d A y calculated using AV18 (dashed 
line). The calculations are compared to data from ref. p4| at E[ ab = 5,7,9,10,12,16,18 



MeV, and from ref. [£5[ at = 22.7 MeV. As expected, Coulomb effects are appreciable 
in all the energy range. Below 18 MeV the effects are appreciable at the maximum. Above 
18 MeV the shape of A y changes and a clear minimum appears, where Coulomb effects 
can be observed. The order of magnitude of these effects is 15%. Instead, 3NF effects 
are not so important. This is a characteristic of the Urbana potential that modifies the 
quartet P-waves (which produces the main contribution to A y below 30 MeV) in such a 
way that there is a cancellation among the different contributions, so the global effect on 
the observable is small. A similar analysis holds for iTu, shown in Fig. 5, since these two 



observables have rather similar structures. The calculations for %T\\ are compared to data 
from ref. @ at E lab = 5, 7, 9 MeV, from ref. JEJ at E lab = 10, 12, 16.5, 22.7 MeV and from 



ref. at = 28 MeV, taking care again that at 16 MeV the comparison is to data 
obtained at a slightly different energy (16.5 MeV). As a difference between the proton and 
deuteron analyzing powers, we observed that Coulomb and 3NF effects are of the same size 
at the minimum of %T\\ above 16 MeV. 

In Figs. 6-8 the tensor observables T 20 , T 2 \, T 22 are given, respectively. As before, the three 
curves correspond to the p-d calculated using AV18 (solid line) and AV18+URIX (dotted 
line), and the n-d T^- calculated using AV18 (dashed line). The calculations are compared 



to data from ref. at E lah = 5, 7, 9 MeV, from ref. g§ at E lab = 10, 12, 16.5, 22.7 MeV 



and from ref. |22| at Ei ab = 28 MeV. As a general trend, the agreement with the data for the 
tensor observables is better than for the vector observables. Coulomb effects are appreciable 
in the three observables at low energies. As the energy increases the inclusion of the Coulomb 
interaction in the analysis of the tensor observables is less important, mostly for T 20 and 
T 22 . The case of T 21 is of particular interest since Coulomb effects are still appreciable at 
28 MeV. This observation suggests that comparisons of d-p data to calculations where the 
Coulomb interaction has been neglected should be done with caution. The effect of the 3NF 
is somehow contradictory since in some cases its inclusion improves the description of the 
observables but in other cases it does not. For example, a net improvement is obtained in 
the description of the minimum of T 22 below 12 MeV. Also the maximum of T 20 and T 21 is 
better described with the AV18+UR potential. Conversely the description of the minimum 
of T 20 and T 21 is better described by the AVI 8 potential alone. The case of T 21 is again 
of interest since 3NF effects seem to be bigger in this tensor observable than in the others. 
Unfortunately experimental data for T 2 \ are not available at all the energies. 

In order to give a quantitative estimation of the agreement between the theoretical 
calculations and the measurements, the Xn f° r the polarization observables is presented in 
Table III. In the first row of the table, the x% * s gi ven with respect to a recent measurement 



performed at 1 MeV jy] and in the second row with respect to the measurements of ref. 
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at Ei a b = 3 MeV. These two energies are below the DBT and are useful for analyzing the 
trend of \ 2 starting at low energies. By inspection of the table, the manifestation of the A y 
puzzle is evident since the Xn f° r the vector observables is a few hundreds at low energy. 
Above 18 MeV A y and iTn change shape being closer to the shape of T21 with a pronounced 
minimum followed by a maximum. After that energy the values of Xn decrease in such a 
way that, at the last energy, vector and tensor observables have similar values which are of 
the order of one tenth. These final values are comparable to those ones obtained recently in 



ref. [441 at 135 MeV. 



VI. CONCLUSIONS 

In the present paper we have studied p-d elastic scattering above the DBT, up to E^ = 
28 MeV. The differential cross section and the total breakup cross section as well as the vector 
and tensor analyzing powers have been calculated using one of the modern NN interactions, 
namely the AV18 potential. In order to evaluate 3NF effects the three-nucleon potential 
of Urbana has been taken into account. The effects of the Coulomb interaction has been 
considered in the framework of the complex Kohn Variational Principle. 

The internal part of the p-d scattering wave function has been expanded in terms of 
the PHH basis. The KVP has been applied to obtain a set of differential equations for the 
hyperradial functions. The set has been solved imposing outgoing boundary conditions at 
a certain value of the hyperradius p = po and expanding the hyperradial functions in the 
region [0, po] i n Laguerre polynomials plus an auxiliary oscillating function. The solution 
should not depend on the value of p providing that for p > p the asymptotic behavior has 



been reached. Such a technique has proved to be adequate since the results from ref. |30| 
has been reproduced as well as the benchmark of ref. |35| . 



The calculations have been extended to all states and parities with J < 19/2, corre- 
sponding to nine energies up to Eiab = 28 MeV. The elastic 5*-matrix has then been used 
to calculate the observables of interest and compare them to the data. Moreover, the corre- 
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sponding observables for n-d scattering, where the Coulomb interaction is absent, have been 
calculated too. From the analysis of the results some conclusions can be drawn about the 
capability of the AV18 and AV18+URIX interactions to reproduce the data. A quantitative 
measure of the agreement achieved by the theory in the description of the data has been 
given through a x 2 analysis. No appreciable 3NF effects have been observed in the total 
breakup cross section and small effects appear in the differential cross section. However, 
results using the AV18+URIX model produce a lower xii value than those where the AV18 
interaction has been used alone. At low energies this is a manifestation of the correct descrip- 
tion of the 3 He binding energy by the AV18+URIX interaction. But as the energy increases, 
there is a different sensitivity to the 3NF. In particular, at the highest energy considered, 
Eiaf, = 28 MeV, the AV18+URIX differential cross section is above the AV18 differential 
cross section, reversing the order observed at lower energies. This situation becomes much 
more evident for example at Ei a f, = 60 MeV [|J. In order to further analyze this behavior 
in Table II the minimum of the AV18 and AV18+URIX cross sections are compared to the 
data. The minimum of the n-d differential cross section calculated with the AV18 potential 
is also given in order to estimate Coulomb effects. The results of the table are useful for 
analyzing what has been called the "Sagara discrepancy" . 



In the case of vector and tensor analyzing powers the Urbana 3NF has little impact 
below 30 MeV. The centrifugal barrier is still strong and there is not too much sensitivity 
to the short range part of the interaction. Conversely there are important Coulomb effects. 
In order to improve the description of the vector analyzing powers new terms could be 
considered in the three-nucleon potential. In such a case both the vector and the tensor 
analyzing powers should improve as well. 

The present picture of the 3N scattering from the theoretical point of view is the following. 
At low energies there is a large underprediction of the vector analyzing powers whereas the 
differential cross section and tensor analyzing powers are well described. Up to 30 MeV we 
see an improvement in A y and iTu indicating that the A y puzzle is a low energy problem. 
On the other hand, a progressive deterioration in the description of the cross section and 
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tensor observables is revealed through a x 2 analysis. For energies above 30 MeV we can refer 
to the very recent work of ref. [|J and we see that this picture remains essentially the same 
up to very high energies (around 135 MeV). Above this energy a number of new conflicts 
appears. 

At present a few realistic local and non-local NN interactions have been determined 
by accurately fitting the two nucleon scattering observables. All of them give rise to the 
A y puzzle in the three-nucleon system. It seems difficult to derive a new NN interaction, 
which still accurately fits the two-nucleon scattering, and correctly describes the N-d A y 
at low energies fSOfl . So, if this is the case, a possible solution to the puzzle should come 
from an improvement of the currently used 3NF models. Indeed, in accordance with chiral 
perturbation theory, these 3NF's include those terms of larger magnitude. On the other 
hand, the A y puzzle can be solved by rather small changes in certain P-wave phase-shifts, 
which can be obtained by adding a small term to the three-nucleon potential ||45|| . However, 
as the energy increases other discrepancies arise in the polarization observables that could 
have different origins. Hence, since the three-nucleon continuum can be calculated at present 
with great accuracy, it is reasonable to expect that the actual 3NF models can be adjusted 
to describe the 3N data. The necessary calculations will require large computing time, but 
is our opinion that this project will certainly help to understand the long unsolved problems 
in low energy nuclear physics. 
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TABLES 

TABLE I. Phase-shift and mixing parameters for different values of the number M of the 
Laguerre polynomials used in the expansion of the hyperradial functions. The s-wave potential of 
Malfliet-Tjon has been considered. 





n-d at E n = 14.1 MeV 


M 




2 r]o 


% 




A 

4 


104.44 


0.4672 


68.993 


0.9669 


« 

o 


105.33 


0.4663 


68.963 


0.9774 


1 9 


105.42 


0.4658 


68.951 


0.9781 


16 


105.49 


0.4646 


68.952 


0.9782 


20 


105.48 


0.4648 


68.952 


0.9782 


24 


105.48 


0.4649 


68.952 


0.9782 


28 


105.48 


0.4649 


68.952 


0.9782 


ref.[27l 

L J 


105.50 


0.4649 


68.95 


0.9782 


Los Alamos 


105.48 


0.4648 


68.95 


0.9782 


Bochum 


105.50 


0.4649 


68.96 


0.9782 




p-d at E p = 14.1 MeV 


4 


107.37 


0.5006 


71.665 


0.9654 


8 


108.34 


0.4984 


72.615 


0.9799 


12 


108.42 


0.4988 


72.602 


0.9794 


16 


108.45 


0.4984 


72.602 


0.9795 


20 


108.43 


0.4984 


72.604 


0.9795 


24 


108.44 


0.4984 


72.604 


0.9795 
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28 


108.44 


0.4984 


72.604 


0.9795 


ref.[27] 


108.43 


0.4984 


72.60 


0.9795 




n-d at E n = 42.0 MeV 


M 




2 r]o 


% 


4 % 


4 


42.198 


0.4575 


38.218 


0.8917 


8 


41.818 


0.4934 


37.680 


0.9028 


12 


41.147 


0.5009 


37.607 


0.9016 


16 


41.271 


0.5010 


37.724 


0.9027 


20 


41.332 


0.5020 


37.723 


0.9031 


24 


41.340 


0.5022 


37.722 


0.9033 


28 


41.341 


0.5022 


37.722 


0.9033 


ref.[27] 


41.33 


0.5026 


37.71 


0.9034 


Los Alamos 


41.34 


0.5024 


37.71 


0.9035 


Bochum 


41.37 


0.5022 


37.71 


0.9033 




p-d at E n = 42.0 MeV 


M 




2 r]o 


% 


% 


4 


38.048 


0.3988 


39.937 


0.8738 


8 


44.701 


0.4961 


40.554 


0.9103 


12 


43.479 


0.5079 


39.844 


0.9011 


16 


43.618 


0.5059 


39.934 


0.9036 


20 


43.660 


0.5055 


39.947 


0.9043 


24 


43.667 


0.5056 


39.947 


0.9046 


28 


43.667 


0.5056 


39.947 


0.9046 


ref.[27] 


43.65 


0.5058 


39.94 


0.9047 
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TABLE II. The minimum of the n-d and p-d differential cross sections (in mb/sr) at different 
energies, calculated using the AV18 and AV18+UR potential models. Experimental data for the 



p-d cross section are from refs. [22,24,25,4 



Energy 


AV18(n-d) 


AV18(p-d) 


AV18+UR(p-d) 


exp. 


1 MeV 


148.9 


177.5 


170.7 


170.2±1.3 


3 MeV 


92.7 


96.0 


92.3 


91.1±0.7 


5 MeV 


53.1 


56.2 


53.8 


52.7±0.4 


7 MeV 


32.9 


35.3 


33.9 


32.9±0.2 


9 MeV 


21.3 


23.1 


22.1 


21.8±0.2 


10 MeV 


17.3 


18.9 


18.2 


18.0±0.2 


12 MeV 


11.8 


12.8 


12.5 


12.2±0.1 


16 MeV 


6.0 


6.5 


6.4 


6.2±0.1 


18 MeV 


4.5 


4.9 


4.9 


4.7±0.1 


22.7 MeV 


2.7 


2.8 


2.9 


2.89±0.03 


28 MeV 


1.8 


1.8 


1.9 


2.19±0.02 
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TABLE III. x 2 P er datum obtained in the description of the vector and tensor analyzing powers 
at several energies using the AV18 and AV18+UR potentials 



Energy 


potential 


A 

Ay 


iT u 


T20 


T21 


T22 


1 MeV 


AV 18 


2(b 


1 1 o 

liz 


3.5 


A £ 

4.5 


2.8 




A t t -i o i T TT~» 

AV18+UR 


190 


bl 


1.0 


2.5 


0.7 


6 MeV 


AV 18 


313 


nnr 

2U5 


4.8 


b. 1 


1 

IZ 




AV18+UR 


271 


144 


5.4 


11 


2.4 


5 MeV 


AV18 


OI 1 

211 


99 


6.8 


1 
12 


7.8 




A T 7" -| o i T TT"» 

AV18+UR 


186 


59 


26 


3b 


1.5 


7 MeV 


AV18 


303 


on 

90 


1 n 

19 


38 


1 n 

1.9 




4 t j -i o i T TT~» 

AV18+UR 


239 


5b 


40 


81 


4.2 


9 MeV 


AV18 


ono 

292 


lb5 


A O 

42 


70 


38 




AV18+UR 


218 


134 


63 


8b 


7.2 


11) MeV 


AV18 


o o o 

288 


29 


1 n 
10 


b.2 


24 




AV18+UR 


224 


23 


13 


b.l 


7.b 


1 O A It ~\ T 

12 MeV 


AV18 


313 


50 


1 n 

19 




on 

39 




A T T -1 l~> • T TT1 

AV18+UR 


227 


34 


lb 




22 


lb MeV 


AV18 


29b 


on 

80 


114 




70 




AV18+UR 


246 


bl 


139 




48 


1 O A T 

18 MeV 


AV18 


ono 

293 












AV18+UR 


250 










22.7 MeV 


AV18 


78 


89 


44 




24 




AV18+UR 


72 


bl 


59 




17 


28 MeV 


AV18 




19 


10 


7.1 


11 




AV18+UR 




13 


10 


11 


8.5 



29 



Figure Captions 

Fig.l. Phase-shift and mixing parameters (in degrees) as a function of energy, (a) 2 S 1 / 2 
and 4 S , 3/2, their real and imaginary parts are indicated by (0> +) an d x), respectively, 
(b) 2 Pi/2 and 2 -P3/ 2 , their real and imaginary parts are indicated by (0;+) and ( n , x ), 
respectively, (c) 2 D 3 / 2 and 2 D 5 / 2 , their real and imaginary parts are indicated by (O, +) 
and (□, x), respectively, (d) 4 -Pi/ 2 , A Pz/2 and 4 -P5/2, their real and imaginary parts are 
indicated by (0>+)> ( D > x ) and (O, *), respectively, (e) 4 -Di/ 2 , 4 -D3/2, 4 -Ds/2 and 4 -D 7 / 2 , 
their real parts are indicated by 0> ^ and A, respectively. Only the imaginary part 
of 4 -Di/2 is given (+). (f) Mixing parameters f?i/ 2 +, ^3/2+, £1/2- and e 3 / 2 _, their real parts 
are indicated by 0> n > A and y, respectively. The imaginary parts of ei/ 2 _ and e 3 / 2 _ are 
indicated by + and x, respectively. The imaginary parts of 771/2+ and r] 3 / 2+ are close to zero 
and are not shown. 

Fig. 2. The p-d total breakup cross section a b below 30 MeV calculated using the AV18 
interaction. The experimental results of Gibbons and Macklin (open triangles) and 
Carlson et al. (open circles) are given for the sake of comparison. 

Fig. 3. N-d differential cross section up to 28 MeV. Calculations are shown for p-d scatter- 
ing using the AV18 (solid line) and AV18+UR (dotted line) potentials and for n-d scattering 
using the AV18 potential (dashed line). Data are from ref. |24|] at Ei ab = 5, 7, 9, 10, 12, 16, 18 
MeV, from ref. |25[ at E lab = 22.7 MeV and from ref. g| at E lab = 28 MeV. 

Fig. 4. Nucleon vector analyzing power A y up to 28 MeV. Calculations are shown for 
p-d scattering using the AV18 (solid line) and AV18+UR (dotted line) potentials and for 
n-d scattering using the AV18 potential (dashed line). Data are from ref. @] at Ei ab = 
5, 7, 9, 10, 12, 16, 18 MeV, and from ref. at E lab = 22.7 MeV. 



Fig. 5. Deuteron vector analyzing power iTn up to 28 MeV. Calculations are shown for 
p-d scattering using the AV18 (solid line) and AV18+UR (dotted line) potentials and for n-d 
scattering using the AV18 potential (dashed line). Data are from ref. at Ei ab = 5,7,9 
MeV, from ref. |2BJ at E lab = 10, 12, 16.5, 22.7 MeV and from ref. at E lab = 28 MeV. 



Fig. 6. Tensor analyzing power T 2 q up to 28 MeV. Calculations are shown for p-d scatter- 

30 



ing using the AV18 (solid line) and AV18+UR (dotted line) potentials and for n-d scattering 
using the AV18 potential (dashed line). Data are from ref. |M] at E\ ab = 5, 7, 9 MeV, from 
ref. pi at E lab = 10, 12, 16.5, 22.7 MeV and from ref. |22[ at E lab = 28 MeV. 



Fig. 7. Tensor analyzing power T21 up to 28 MeV. Calculations are shown for p-d scatter- 
ing using the AV18 (solid line) and AV18+UR (dotted line) potentials and for n-d scattering 
using the AV18 potential (dashed line). Data are from ref. p4| at E\ ab = 5, 7, 9 MeV, from 
ref. p| at E lab = 10 MeV and from ref. pi at E lab = 28 MeV. 



Fig. 8. Tensor analyzing power T22 up to 28 MeV. Calculations are shown for p-d scatter- 
ing using the AV18 (solid line) and AV18+UR (dotted line) potentials and for n-d scattering 
using the AV18 potential (dashed line). Data are from ref. [^4| at Ei ab = 5, 7, 9 MeV, from 



ref. [251 at E lab = 10, 12, 16.5, 22.7 MeV and from ref. 221 at E lab = 28 MeV. 



31 



180 
135 
90 



cn 

~0 4b 



CO 



+ 

C\J 





-45 
-90 



en 

CD 

Q 




4 



cn 

CD 
T3 



-4 



-8 



-12 



O 



O 



o 



o 



□ 



(a) 

o 



o 

□ □ 



□ 



EL 



(C) 



* 



* 



(e) 



o 

□ 



*^<>^XX XX X X - 



A 



10 20 
E lab[ MeV ] 



O 

□ 



A 



30 



20 



10 



en 

CD 
T3 



(b) 



x □ 



□ 

□ o' 
o 



□ o 



O 



,o 



O 



o 



10 20 
E lab[ MeV ] 



x 



O 




30 



200 




E lab [MeV] 



